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Abstract 



The late stages of the nucleation have been described analytically. 
The approximate solution of the Zel'dowich-Folmer-Frenkel equation 
has been constructed. 



Introduction 

The process of the evolution after the end of nucleation (the main stage) 
and the consumption of the main quantity of the surplus metastable phase 
requires a separate description. This description can be based on the real 
form of the size spectrum [5] or can be done by the asymptotic analysis pQ , 

0- 

The regime will be chosen as the free molecular one. This is done because 
the formation of the exhausted zones around the droplets strongly influences 
the kinetics and the description can not be based on the averaged character- 
istics. The only situation where the diffusion regime can be combined with 
the averaged characteristics is when the embryos go from one cave to another 
and change many exhausted zones. But this situation can be reduced to the 
effective absence of the cave and, thus, the regime will be the free-molecular 
one. It is important because earlier there was the restriction that the size is 
less than the length of the free motion of a molecule divided by the coefficient 
of condensation. Now this restriction is not essential. 

The first stages of the late evolution were investigated in [5]. 
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1 Balance equation 

Define n as the number of embryos with v molecules. The evolution in time 
t is described by equation 

Oti 

— = W + (u - l)n(u - 1) - W-(v)n(v) - W + {v)n(v) + W~(u + l)n(u + 1) 

Here W + is the direct kinetic coefficient, W~ is the inverse kinetic coefficient. 
In the Fokker-Planck approximation the will be a diffusion and the regular 
growth. 

In approximation of the pure diffusion the balance equation will be 

Oti 

— = W + (u)[n(u - 1) - n(v) - n{v) + n(u + 1)] 

and 

— = W + (u) 



dt dv 2 

2 Some estimates 

According to [5] after the regular relaxation there will be a diffusion errosion 
of the spectrum. For this stage one can give some estimates. 

2.1 Diffusion 

One has to use s-scale. Here s is the surface of the embryo 

This is because the kinetic coefficient of absorption W + is proportional to s 

W + = W°s 
Then in diffusion equation will be 

dt du 2 y } 



In the leading term 



d 2 n(v) _ 2 i z d 2 n{v 



v 



dv 2 ds 2 
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Then 



dn(y) 



= W 1 



o 



2 n(u) 
ds 2 



dt 



The same will be for n(s). 
The diffusion scale will be 



Sdiff - t 



1/2 



For the linear size of an embryo 



P = 



.1/3 



we have 



Pdiff 



f l/4 



2.2 Regular growth 

The coordinate of the regular growth will be determined from 

dp = ^ 
dt t t 

Here t t is the characteristic time, ( is the supersaturation defined by 

C = ni/noo - 1 

where n\ is the molecules number density of vapor, is the molecules 
number density of the saturated vapor. 



Then 



Preg 




The behavior of ( is given by 





3 Pc lpc=pMff 
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Then 

2a 

and 

c ~ t-^ 

Then 

Preg ~ t 3 ' 4 

We see that 

Pdiff ^ Preg 

and it is clear that the main quantity of substance will be in the tail growing 
regularly. 

Approximately we assume that at the period when the tail will be essential 
the following approximation W° ~ const takes place. 



3 The form of the tail 

We determine the form of the tail. Imagine that at big sizes there is a 
dominating regular growth, at small sizes of the critical order there is a 
dominating diffusion. The boundary will be marked by the index b. 

Certainly, at the very small sizes there is a regular dissolution. But this 
region is simply negligible and can be expelled by the transition of a zero 
point. 

Then for p > Pb 



h> t t 

Then the distribution function / is given by 

f(p,t) = fb(Pb,t') 

To get fb one can use the diffusion approximation according to [7], [5] 

, A ( (s b -s c ) 2 ^ 

In the absolutely stationary approximation one can put the normalizing con- 
stant 

A = l 
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Due to the dissolution this constant will be another and this question will be 
solved later. 

Now we consider 

Ofe - s c ) 2 



4W°t' 

As an estimate one can take 

Pb~ Pc = 2p c 

Then 



AW°t' 4W°t' 

If we take p c ~ t 1 / 4 (more precise it is necessary to take p c ~ W°t 1 / 4 ) 
then 

(s b - s c ) 2 ^ 64 
4W°t' ~ 4 

Then 

/& ~ / ° exp(— 16) 
v4l¥ 

4 More precise formula 

One can use another more precise equation (see [7], [5]) 

, A> , (s~s c0 ) 2 A (s + s c0 ) 2 

fb = rTTTTTT eX Pl TrxrTU ) 77T777P eX P( 



Here s c o marks the initial position of a critical surface. 
After the transformations we get 

, A s 2 , s 2 u , 2ss c0 . , 2ss c0 . , 

; = v*pf exp( -W exp( -W [exp( W - exp( -W ] 

One can see that 

ss c0 /AW°t < 1 

Really, since 

and s co is fixed we come to 

ss c0 /AW°t ~ r 1/2 < 1 
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Then after the decompositions of the corresponding exponents one can come 

tO 2 2 

So, 

/ ~ yexp(-y 2 ) 

with 

s 

y 



V4WH 

The maximum of this distribution is attained at 

y m = V2/2 

This value will be close to the critical value 

Vm ~ Vc 

For 

Vb = ay c 

with some parameter a we get 
with constant f3. 

5 Dissolution of the spectrum head 

The dissolution of the head determines the form of the tail. 
At first there is a pure dissolution of a gaussian and 

p c ~ const 

Then 

f b ~ -jJ=eM-(s b - s c ) 2 /(4W°t)) 

or 

h ~ ^L=exp(-9^V^o/(4^ t)) 
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For the behavior of f b we get 

with a positive parameter 7 = const. 

The last function is positive, it has an evident asymptotic behavior /& — > 
t -1 / 2 at big i. At small t it goes to zero. Thus, it has a maximum. 

To see the formation of the tail one can approximately write 

P~ PcO = —{t-t) 

H 

where Co is some initial supersaturation. Here we assume Co to be a constant 
which is a very strong assumption. So, below all results will have the power 
of stimates. 

The maximal size pu m is 

Plim = PcO + Co^Ai 

Then having put the starting moment at zero one can get 

t' = (p- pl im )t t /( 
The distribution at the tail looks like 

/p = /M = (|4 1/2 exp(-^) 

as a function of 

Z — Plim — P 

Now we shall calculate the quantity of substance in the tail G ta ii- We get 
Gu#= f^A-r-^ -) 1 / 2 exp(- 7 — ) dp 

J3p c t t {pum - P) UyPlim - P) 

Having introduced 

Plim ~~ P > / 1 

y = — - — e = iCo/U 

we get 

Plim- 3 Pc j 



One can see that G tai \ grows very fast. The asymptotics is 

Gtail -> Pilm ~ 

But this behavior will take place only until 

ss c ~ 4vy°t 

After this moment of time the boundary condition will be another 

f\ Pb = A exp(-const) 

or at larger scales of time 

Alternative approach 

At the further stage 

So, we see that the variation is very small. The drift of p is rather small 
also 

p ~ J (dt ~ t 3/4 

because of 

The subintegral function ^ tai j in the quantity of substance 

Gtail = J dtaildp 

in the tail grows like 

gtau ~ t 9/4 

It grows fast. 
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6 Impulse regime 



Since G ta ii grows very fast one gas to analyze this behavior. 
At first 

Gtail G total (2) 

where G to tai is the total quantity of substance in droplets. 

The main role here is played by diffusion and the solution can be described 
with the help of combinations of Gaussians. This model will be called as the 
diffusion-regular model. It is investigated in [5]. 

Since g ta n grows fast it means that the monodisperse approximation for 
Gtaii is quite suitable 

Gtail ~ PeffNeff 

Here N e ff is the effective number of droplets in the tail and p e // is their 
coordinate. Certainly, N e ff is not a constant value, it grows. 
At t determined by 

Gtail ~ Gtotai — Gtail 

the tail begins to eat the head of the spectrum. Here the tail is non-essentially 
decreasing one. It occurs very fast and the spectrum eats at first the head 
and then it begins to eat the the rather flat beginning of the tail. 

This process can be easy described by the regular growth of the tail and 
by the regular dissolution of the eaten part of the spectrum. 

Then only the end of the tail will exist. This end can be then treated as 
a monodisperse spectrum. Later this part will be dissolved by diffusion and 
the process repeats. 

We shall call this process as the "impulse condensation". 

This part is described in [5]. 

7 Short tail 

A question whether the process of impulse condensation repeats arises here. 
From the first point of view it seems that this process will infinitely repeat. 
But there is one objection. 

The function f(t) = exp(— r y/t)/y/i described above has a maximum. 
This maximum takes place &tt/ r y = y = y m ~2 and for y essentially less 
than y m there is a rapid decrease to zero. 

Certainly, after every cycle the time t has to be shifted. 
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The crucial point is whether at y m the condition (j2J) is observed. An 
alternative case 

Gtail ~ G total — Gtail 

at y m leads to another further evolution. Here the tail is essentially decreas- 
ing. This case will be called as "the case of adjusted tail". 

In this section one can study the evolution on the base of the regular 
growth. This radically simplifies the situation. 

Here we shall give the main ideas of description in this situation. The 
vapor consumption by the head of the tail can eat only the same head of the 
tail. There is no other parts to eat. Then we come to the following balanced 
picture: 

• The task is to determine p& = 2 + 3p c , p c = 2a/31n(£ + 1) « 2a/3( 
where a is a renormalized surface tension. 

• The size distribution for p > p& is known. We have 

f(p,t) = f(x) 



for 

z = pn m (t). Then 
For f(x) we have 



x = z — p 



r*-Pb „ 
Gtaii = - x) f(x)dx 



f(x) ~ exp(--)/Vx 

X 



The balance equation is 

Gtaii + G head + noo(C + 1) = const 

where Ghead is the quantity of the substance in a region 0.7p c < p < pb- 
In the region p < 0.7 p c one can simply neglect the substance. 



In the region 0.7 p c < p < p one can write the purcl diffusion equation 
(JTJ) with boundary conditions 

f{p,t) = 



-^.e. without the regular growth and with the constant coefficient of diffusion. 
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at p = 0.7 p c 

f(P,t) = ftail(x) 

at p = p b 

• In the region p < 0.7 p c one can simply neglect the evolution. 

8 Further simplification 

Already the last system of equations can be solved when we know f(x). The 
function f(x) has to be known from the initial conditions. These conditions 
are: 

• The direct result of the nucleation period [I]. 

• The result of the impulse regime [5]. 

To show the result in a most simple way one can use some simplifications. 
Further simplification is the following 

• To write the balance equation as 

Gtail = fixed 

• To use that 

Gtail ~ Ntailpb 

where N ta u is the effective (proportional) number of the droplets in the 
tail. 

The function f ta u is one of the possible realizations, which has a certain 
disadvantage. Really, for p > pu m the distribution is zero. Meanwhile, the 
diffusion in the region p > pb inevitably leads to the big droplets. Assume 
that the behavior of the gaussian at big arguments can be approximated as 

/ ~ exp(-7'p) 

with some parameter 7'. 

The same asymptotic can be derived from f ta u by the steepest descent 
method. 
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So, the last approximation is very fruitful. It allows to calculate the 
integrals for G ta ii analytically. 

One can make some new simplification. Since the spectrum is rather 
sharp, it means that 

dp b dp c 
dt dt 



which leads asymptotically to 



-1/2 



Since the diffusion width is Sd ~ t 1 ^ 2 , pd ~ t 1 / 4 we see that the diffusion 
does not radically change the character of the process. However, since the 
Gaussian does not have a finite support the droplets with the greatest sizes 
which will be the main in further consumption appears due to diffusion. This 
is one more point for the approximation /. 

The constant 7 is not a precise constant but a slowly varying function. 
This effect can be taken into account by standard methods. 

An approximation / is so simple that we have no need to use the monodis- 
perse approximation. One can calculate G ta u directly by integration 



poo 

Gtaii ~ / {z - x) 3 exp(->y'x)dx 



and the integral can be taken analytically. 



9 Initial asymptotic solution 

Since we came to the balance equations typical for the asymptotic solutions 
of Lifshic and Slyozov [1] it is necessary to analyze this approach. 

As a regular law of droplets growth one has to take the precise expression 

dp 6cr 1 1 

dt 9t t p c p 

where a is the renormalized surface tension. 
For 

u = p/p c 
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one gets 

du 6(7 , 1 1 . u dp c 
dt 9t t p c p c p p c dt 
In the reduced coordinates 

^ = (1 - l/u) -7« 
dpc 

Then to escape the violence of the substance balance one has to observe 

dt ~ 

for all u. 

The last relation will take place for every time r increasing in time since 

du du dt 
dr dt dr 

Namely 

r = 3\n(p c /p c (t = 0)) 

was chosen by Lifshic, Slyozov (LS). 

The last requirement is very strong. In reality, one can imagine the situ- 
ation where this requirement is not valid. In the general case it is necessary 
that this requirement takes place only in the integral sense. But in LS theory 
this requirement has to be valid at every moment. 

LS required that 

( du \ - n 
max {u} {—) = r -> 

Meanwhile the real condition should be the following: 

• If r has a power asymptotics (asymptotics as an argument in some 
power) then this asymptotics has to be zero 

But the situation when r has no asymptotics is also quite possible. For 
example, the case of oscillations drops out of attention here. 

The crucial supposition of LS theory is the existence of the asymptotics 
for r (in the class of constants). 
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The recipe of LS is to consider r as the monotonous function and then it 
is reasonable to put r = asymptotically. 

It is necessary to stress that the condition du/dt < for every u leads to 

dp p dp 



for every p or 



- < 
at p c dt 



dlnp ^ d\np c 



dt dt 

for every p. This is very important relation. 

From the last inequality it follows that every p can not escape and will 
inevitably be dissolved. But at some u there will be no dissolution. Let it 
be uq. Requirement r = leads to 

u = 2 

9tt _ A 



n 2 dpc 



The last equation can be easily integrated which leads to 

- Ap~ l = const + t^- (3) 

but one has to note that if we observe the validity of the previous equation 
only in the averaged sense we get approximately the same integral law. So, 
one can not regard the law of 4/9 (this is the characteristic coefficient in the 
diffusion case, here the regime is the free-molecular one) as an experimental 
justification of LS approach. 

The law ([3D can be easily justified. Really, if there will be another power 
asymptotic then immediately we come to the violence of the balance condi- 
tion. 

One has also to stress that it is absolutely impossible to differentiate (j3J), 
the result will have nothing in common with a real situation. This fact is 
typical for such asymptotics. 



10 The form of the spectrum in LS theory 

It is clear that the evolution of the system is governed by the spectrum of the 
droplets sizes. This function determines the behavior of the supersaturation 
in the system. 



14 



In frames of LS approach the behavior of p c is already known. Then it 
is possible to solve the inverse problem: to reconstruct the spectrum on the 
base of the behavior of supersaturation. The sense of this way of derivation 
is wrong, but technically it is quite possible. 

Under the regular law of growth one can find for the distribution function 

0(t,m) = -£(r - t(u))/v u 

where 

v u = du/dr 

is already known and 

ru 

t{u) = \ dujv u 
Jo 

is also known and £ is some arbitrary function which is going to be deter- 
mined. The balance equation allows to determine the form of £ which solves 
the problem. 

To solve it in the simple form r is chosen as 

r = 31n(p c (t)/p c (t = 0)) 

. Then the balance equation in a closed system will be 

exp(r) / M 3 0(r, u)du = 1 
Jo 

which can take place only if 

£(...) = exp(...) 

The spectrum is zero for u > -u and continuously goes to zero for u — > 
Mo — 0. So, the essential part of the spectrum ill be dissolved at finite times 
r and t. Here it is necessary to stress that now it is clear that the deviation 
of r from zero at some finite interval (it is quite possible) can lead to the 
absolutely another form of the spectrum. So, the form of the spectrum is 
determined only with the help of the strong supposition of LS theory about 
the behavior of the critical size, i.e. about the supersaturation. 
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11 Correction for 7 



The limit value r = is not acceptable even in LS theory. The next step 
leads to correction of r. Ordinary, the correction is many times greater then 
the precedent value of r which allows Osipov and Kukushkin [2] to speak 
about the non-uniform character of LS decompositions. 

In frames of LS theory it is possible to write a correction term and to 
fulfill the next step of calculations. Osipov and Kukushkin (OK) [2] follow 
another way. Certainly, since the result in the first approximation is known 
it is possible to rewrite the formulas in such a way that 7 will be zero and 
the first approximation has to inevitably include the correction term in LS 
approach. Since the conclusions of OK theory are important it is necessary 
to consider it. We omit all details which can be founco in [2]. 

Having introduced 

1 [* 

and redefined u = p/2p c one can come to 

du (u — l) 2 

+ 7« 



dr u 

Now the condition r = corresponds to 7 = and every correction for 7 
will be giant in comparison with 70 = 0. Then the asymptotics here has to 
include the correction term already in the main order. 
The approach with correction gives 



7 = 1/(4 



_2^ 



T 



The substitution 

v — (1 — u)~ l — In \u — 1| 
brings the rate of growth to 

^ = i + 7(i + V0 2 

where ip satisfies 

%p + In |^| = v 



2 The initial total number of embryos is found with an error from [3] because the asymp- 
totic analysis [3] is not valid, it is necessary to use [1]. 
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and asymptotically 

dv 1 o 

Tt =1 + ^ 

The initial asymptotics in LS theory corresponds to 7 = and then 
asymptotically 

v — t + const 

The initial asymptotics in OK theory corresponds to 7 = l/(4r 2 ). Then 
asymptotically 

v -> 2r 



12 Correction for the distributions 

The distribution function can be in frames of LS and OK theories found from 

f{u,T) = ac{u^))— 

where C is the integral of the law of growth. 

Concrete calculations will give the following results 

• In LS theory 

v — t — const 

or 

exp(t> — r) = const 

i.e. 

Cls = z exp(-- r) 

1 — u 1 — u 

is the integral of evolution Then £ is 

i~C 2 LS at C LS >0, £-0 at C LS <0 
Then the spectrum is 

2u . 2u , 
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In OK theory 

v — 2r = const 

or 

exp(t> — 2t) = const 

i.e. 

CW = t exp(-- 2r) 

1 — u 1 — it 

is the integral of evolution Then £ is 

£ ~ CW CW > 0, £ ~ at C k < 
Then the spectrum is 

0Oft- = 7^ ^exp( 



and it essentially differs from LS theory. 

13 Further corrections 

The essential difference of spectrums is the striking feature of the LS theory. 
But it is rather easy to see that such a feature will be natural for all further 
corrections. 

Suppose that the integral of evolution at the previous step is established. 
Then we can reformulate the rate of growth as 

du 

- = F{u) - 1U 

where F is a known function. 

Here we use for simplicity t instead of r. 
In the zero approximation (LS theory) 

F ~ — u ^ + 1/2 

in the first approximation^] (OK theory) 

F= {u-lf/u 



3 With a renormalization 
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etc. 

Let w be the solution of the equation with 7 = 0, i.e. 




— t + const 



Then at the previous step 
At the current step 

dw 
~dl 

So, w is a straight analog of v. 
There can be two situations: 

1. Asymptotically w — > t. 

2. Asymptotically w — > l(t) ^ t. 

In the first situation the spectrum remains absolutely the previous one. 
There is absolutely no corrections. 

In the second situation the correction is essential. Instead of w — t as the 
integral of evolution one has to use w — l(t). It means that instead of w one 
has to use l^~^(w) where is the inverse function. 

The difference in spectrums is striking. 

So, the alternative is to have no corrections or the striking corrections. 
This is the consequence of the use of Lifshic-Slezov variables and a certain 
disadvantage of LS and OK approaches. 

14 Application of asymptotics 

The ideology of LS and OK asymptotic analysis is one and the same: 

• The asymptotic of some given function of the supersaturation in the 
class of powers is prescribed. 

• The balance condition at asymptotics leads to the determination of the 
supersaturation as the function of time. 



dw 
~dt ~ 

= 1 + ju(w) 
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• The form of the size spectrum is reconstructed on the base of the given 
supersaturation. 

Meanwhile, under the regular growth of the embryos the natural sequence 
of actions is the following 

• The size spectrum is given from initial conditions 

• The dissolution of the size spectrum together with the balance equation 
determines the behavior of the supersaturation 

• The tail of the spectrum at big sizes determines the asymptotics of the 
supersaturation and of coefficients in the regular law of the droplets 
motion. 

We see that from the last point of view the LS and OK theories are 
inconsistent. Then why they correspond to the asymptotic behavior found 
from experiment? 

At first one has to mention that the accuracy of experimental measure- 
ments is not high: even the striking difference between LS and OK drops out 
of experiments. Hence, we have to conclude that the experiment gives only 
the approximate form of the size spectrum with maximum and two different 
rather short wings. 

Now we shall see that the form of the size spectrum obtained in LS and 
OK theories is rather typical at least approximately. 

The tail of the Gaussian can well approximated by 

£ as = exp(- const * p) 

with some constant const. We shall call this constant A and suppose it to 
be a slow function of p to have 

A(p + A- 1 )-A(p)«A(p) 

Then the exponential approximation will be valid at least at the essential 
part of the size spectrum. 

These transformations are absolutely adequate to the standard approach 
of the steepest descent method. 

The same approximation can be established in tt-scale. The constant can 
simply cancelled by renormalization. Then 

£ as ~ exp(-tt) 
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Now we shall approximate t(u) for two values U\ and m 2 at the essential 
part of the size spectrum. The simplest approximation is the following 

t{u x ) - t(u 2 ) = < —t— ( u i ~ u 2) 

which is the linear connection. Then the exponential form of the spectrum 
over t{u) will be conserved 

£ as ~ exp(-r(w)) 
Since £ has to be the function of r — r{u) we come to 

f = exp(r - t{u)) 

which lies in the base of LS and OK reconstruction of the size spectrum. 

To get the spectrum one has to divide £ on v u . The last value is deter- 
mined by the regular law of growth. Certainly, one can not guarantee that 
the parameter 7 corresponds to r = 0. When it really corresponds to r = 
there will be LS asymptotic. When 7 corresponds to some small r less than 
zero then there will be an extremum associated with the minimum of v u . 

One can see that since r is rather small in the absolute value there is a 
maximum of near uq. One has to note that in the experiment there is no 
evident way to determine the critical size. Then there is no way to determine 
the position of the size spectrum, but only it's form. 

Since A is not a true constant one can say about the quasistationary 
solution with a smooth variation of parameters and a form of the spectrum. 

These considerations evoke the ideas associated with the square approxi- 
mation of the rate of growth du / dr near the maximum which were presented 
in ®. 



15 Direct determination of 7. 

Although the precise asymptotic can not be directly used to determine the 
evolution of the system one can suggest another approach. 

Now we return to the description of the evolution. 

The rate of growth can be written as 

du - , 
— = F(uMr)) 
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The integration of this equation gives f(u). Symbolically this can be written 

/ 



as 

du 



where 

du 

v u = ~r 

dr 

The problem of integration of the rate of growth is very important and ana- 
lytically it is difficult to do, but here we suppose that it is fulfilled. 
The balance equation has to written as 

Gtotai = const 

It is possible to neglect the surplace mother phase. It can be reduced to 

exp(r) / u (f)(u,T)du — const 
Jo 

where <fi is the distribution function, uu m is the maximal size. 
The regular law of growth leads to 

t(r-T{u)) 

<f)(u,T) = 

-v u 

Then the balance equation will be 

r u um u 3 £(t - f(u)) , 

exp(r) / du = const 

Jo —v u 

But the last equation is not an equation on the function £ as in the LS 
and OK theories but an equation on 7. The function £ is known from the 
initial conditions for this equation. This is the crucial difference between our 
theory and LS or OK theory. The dependence 7(7") stands in v u and in f(u). 
Moreover in f(u) it is at least the dependence on 7 in all preceding moments 
of time. So, the last equation is very complex and one gas to suggest some 
methods of the approximate solution. 



16 Approximate solution 

When 7 is really smaller than 70 corresponding to r = it is possible to 
see that some approximations for 7 produces approximations for evolution. 
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When the main influence is ensured by the supercritical tail then the condi- 
tion 

Gtaii = const 

is rather productive. 

One can approximately say that in the region p < pb the rate of growth 
is very small. 

Another possibility is to use for du/dr an approximation which can be 
integrated analytically. Namely du/dr has to be used because for u the 
boundaries of the nearcritical region and supercritical region will be fixed. 
For 

du _ x _ 

- = (1- M )- 1U 

we use the following constructions: 

• We give definitions 

F 1 = (l-u~ 1 ) 
F 2 = r 
F 3 = 1 — ju 

• The first interval will be [u 2 , oo[. The value of u 2 is the root of equation 

F 2 = F 3 

At the first interval one can approximately write 




The last equation can be integrated analytically 
• The second interval is [1*1,1*2] where U\ is the root of equation 

F 2 = Fi 

At the second interval one can approximately write 




The last equation can be integrated analytically 
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• The third interval is [0, u±\. At the third interval one can approximately 
write 

du _ 
Tr~ Fl 

The last equation can be integrated analytically 
Summarizing one can approximately write 

with an analytical solution. Then we know f{u) and can write a closed 
equation on 7. 

This approach will be called as the "regular- regular model". 
In the general situation allowing LS asymptotics 

dujdr = Fx — 7« 

and the procedure remains the same. 

Another possibility is the quasistationary approximation. Here 7 is sup- 
posed to be a local constant and the law of growth is integrated 

r du 

T[U) = / — — 

Then we know t(u) and £(r — t(u)) as function of u,r. Then the balance 
equation is the closed equation on 7. 

Define as r total the time of dissolution from the size u m , where dujdr 
attains maximum, up to zero. It is clear that the necessary condition for the 
applicability of the quasistationary approximation is 

Vtotal(l{j)) - Ttotal(^(r + T tota i (^(t) ))) | ^ 
Ttotal(l(r)) 

The difference between LS and OK theories symbolizes that at the asymp- 
totics r = the quasistationary approximation can be hardly applied at the 
very end of the process. But it can be fruitful at some earlier periods. 
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17 Combination of approaches 

Here the regular-regular model and the diffusion-regular models of solution 
were presented. Now it is necessary to decide what model has to be used. 

The main object will be the support of spectrum. We shall define the 
essential support of the spectrum by the following way: 

• For Gtotai we define the subintegral function g by 

G total = J dpg p 

or 

Gtotai = J dug u 

• We define the maximum of g and 

pmax = &rg(max(g)) 
u max = &rg(max(g)) 

• Then we determine the region of support by 

dug u = 5/2 



/ dug u = 5/2 

J u r 

or by 

rpi 

/ dpg p = 5/2 

J 

POO 

/ dpg p = 5/2 

J Pr 

with some small parameter 5. 

This region can be also determined by conditions 

g{u{) = g(u r ) = 5g(u max ) 



25 



If the essential support belongs to the near-critical region then we have 
to use the diffusion-regular model. In the opposite situation we have to the 
regular-regular model. 

The impulse regime contains the sequential application of the diffusion- 
regular model and the regular-regular model. Here the evolution during 
the period described by the regular-regular model occurs in a very simple 
manner - the essentially supercritical droplets determine the behavior of the 
supersaturation. So, there is no need to use the complete formalism of the 
regular-regular model. 

In the regime of adjusted tail one has to use the regular-regular model 
which can not be simplified. 

To see the times when the regular-regular model will be no longer valid 
it is possible to apply the diffusion operator to the solution obtained in the 
regular-regular model. When the difference between this solution and the 
result of the application of diffusion operator will be essential it means that 
it is necessary to introduce the diffusion corrections. 

The presence of diffusion corrections will be very essential at the tail of 
the size spectrum and namely this tail will govern the evolution later. This 
tail can be hardly described by the diffusion equation because there is no 
statistics - only few droplets of big sizes determine the evolution. 

Ordinary the process of nucleation results in existence of several big 
droplets with stochastic sizes and namely their competition in growth deter- 
mines the latest stage of the process. The sizes of these droplets are rather 
arbitrary. So, the kinetics of interaction will be also unstable and arbitrary. 
This effect forms a matter of a special publication. 

Here it is necessary only to add that the number of droplets strongly 
diminishes in time and to see the impulse regime and later the regime of 
adjusted tail it is necessary to have a nucleating system of huge sizes. 
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